home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
CD ROM Paradise Collection 4
/
CD ROM Paradise Collection 4 1995 Nov.iso
/
program
/
nrpas13.zip
/
RATINT.DEM
< prev
next >
Wrap
Text File
|
1991-04-29
|
746b
|
36 lines
PROGRAM d3r2(input,output);
(* driver for routine RATINT *)
CONST
npt=6;
eps=1.0;
TYPE
glnarray = ARRAY [1..npt] OF real;
VAR
dyy,xx,yexp,yy : real;
i : integer;
x,y : glnarray;
FUNCTION f(x,eps: real): real;
BEGIN
f := x*exp(-x)/(sqr(x-1.0)+eps*eps)
END;
(*$I MODFILE.PAS *)
(*$I RATINT.PAS *)
BEGIN
FOR i := 1 to npt DO BEGIN
x[i] := i*2.0/npt;
y[i] := f(x[i],eps)
END;
writeln('Diagonal rational function interpolation');
writeln;
writeln('x':5,'interp.':13,'accuracy':14,'actual':12);
FOR i := 1 to 10 DO BEGIN
xx := 0.2*i;
ratint(x,y,npt,xx,yy,dyy);
yexp := f(xx,eps);
writeln(xx:6:2,yy:12:6,' ':4,dyy:11,yexp:12:6)
END
END.